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Complex patterns generated by the time evolution of a one-dimensional digitalized coupled map 
lattice are quantitatively analyzed. A method for discerning complexity among the different patterns 
is implemented. The quantitative results indicate two zones in parameter space where the dynamics 
shows the most complex patterns. These zones are located on the two edges of an absorbent region 
where the system displays spatio-temporal intermittency. 
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This century has been told to be the century of complexity ,JJ. Nowadays the question 'what is complexity?' is 
circulating over the scientific crossroads of physics, biology, mathematics and computer science, although under the 
\ present understanding of the world there seems to be no urgency to answer the above question. However, many 
different points of view have been developed to this respect and hence a lot of different answers can be found, at 
present, in the literature |2j. 

It should be kept in mind that in ancient epochs, time, space, mass, velocity, charge, color, etc. were only per- 
ceptions. Only after that they became concepts, different tools and instruments were invented for quantifying those 
perceptions, and, finally, only with numbers the scientific laws emerge. In this sense, if by complexity it is to be 
understood that property present in all systems attached under the epigraph of 'complex systems', this property 
should be reasonably quantified by the different measures that were proposed in the last years. This kind of indica- 
tors is found in those fields where the concept of information is crucial. Thus, the effective measure of complexity 
0, the thermodynamical depth J3|and the simple measure of complexity [j| come from physics and other attempts 
I . such as algorithmic complexity 001 > Lempel-Ziv complexity and e-machine complexity Q arise from the field of 
' computational sciences. 

In particular, taking into account the statistical properties of a system, an indicator called the LMC (LopezRuiz- 
Mancini-Calbet) complexity has been introduced 10]. This magnitude identifies the entropy or information stored in 
a system and its disequilibrium i.e., the distance from its actual state to the probability distribution of equilibrium, 
as the two basic ingredients for calculating its complexity. If H denotes the information stored in the system and D 
is its disequilibrium, the LMC complexity C is given by the formula: 
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C({ Pi }) = H({ Pi }) ■ D({ Pi }) = 

= -k (E^ftiogw) • (Eilt (Pi - i) 2 ) (i) 

where {pi}, with Pi > and i = 1, • • • , N, represents the distribution of the N accessible states to the system, and k 
is a constant. 

As it can be straightforwardly seen, the LMC complexity vanishes both for completely ordered and for completely 
random systems as it is required for the correct asymptotic properties of a such well-behaved measure. Recently, it 
has been successfully used to discern situations regarded as complex in discrete systems out of equilibrium [111 Il2l 
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As an example, the local transition to chaos via intermittency in the logistic map, x n +\ = Xx n (l — x n ) presents 
a sharp transition when C is plotted versus the parameter A in the region around the instability for A ~ At = 3.8284. 
When A < At the system approaches the laminar regime and the bursts become more unpredictable. The complexity 
increases. When the point A = At is reached a drop to zero occurs for the magnitude C . The system is now periodic 
and it has lost its complexity. The dynamical behavior of the system is finally well reflected in the magnitude C (see 
Fig. 3a of Ref. [13). 

When a one-dimensional array of such maps is put together a more complex behavior can be obtained dependin g on 
the coupling among the units. Ergo the phenomenon called spatial-temporal intermittency can emerge |l8lll9ll20ll2l) . 
This dynamical regime corresponds with a situation where each unit is weakly oscillating around a laminar state that 
is aperiodically and strongly perturbed for a traveling burst. In this case, the plot of the one-dimensional lattice 
evolving in time gives rise to complex patterns on the plane. If the coupling among units is modified the system 
can settle down in an absorbing phase where its dynamics is trivial |22l |23j and then homogeneous patterns are 
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obtained. Therefore an abrupt transition to spatio-temporal intermittency can be depicted by the system |2J, |25j 
when modifying the coupling parameter. 

In this letter we are concerned with measuring C in a such transition for a coupled map lattice of logistic type. 
Our system will be a line of sites, i = 1, . . . , L, with periodic boundary conditions. In each site i a local variable xf 
evolves in time (n) according to a discrete logistic equation. The interaction with the nearest neighbors takes place 
via a multiplicative coupling: 

x ^ = (4-3pX?)x?(l-x?). (2) 
The variable X™ is the digitalized local mean field, 
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2 ( x i+l + x ?-l) 
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with nint(.) the integer function rounding its argument to the nearest integer, and p is the parameter of the system 
measuring the strength of the coupling (0 < p < 1). 

There is a biological motivation behind this kind of systems pfiL l27| . It could represent a colony of interacting 
competitive individuals. They evolve randomly when they are independent (p = 0). If some competitive interaction 
(p > 0) among them takes place the local dynamics loses its erratic component and becomes chaotic or periodic in 
time depending on how populated the vicinity is. Hence, for bigger Xf more populated is the neighborhood of the 
individual i and more constrained is its free action. At a first sight, it would seem that some particular values of p 
could stabilize the system. In fact, this is the case. Let us choose a number of individuals for the colony (L — 500 for 
instance), let us initialize it randomly in the range < Xi < 1 and let it evolve until the stationary regime is attained. 
Then the black/white statistics of the system is performed. That is, the state of the variable x t is compared with the 
critical level 0.5 for i = 1, . . . , L: if x.^ > 0.5 the site i is considered white (high density cell) and a counter N w is 
increased by one, or if Xi < 0.5 the site i is considered black (low density cell) and a counter Nb is increased by one. 
This process is executed in the stationary regime for a set of iterations. The black/white statistics is then the rate 
(3 = Nb/N w . If (3 is plotted versus the coupling parameter p the Figure^is obtained. 

The region 0.258 < p < 0.335 where j3 vanishes is remarkable. As stated above, (3 represents the rate between 
the number of black cells and the number of white cells appearing in the two-dimensional digitalized representation 
of the colony evolution. A whole white pattern is obtained for this range of p. The phenomenon of spatio-temporal 
intermittency is displayed by the system in the two borders of this parameter region (Fig. [21 . Bursts of low density 
(black color) travel in an irregular way through the high density regions (white color). In this case two-dimensional 
complex patterns are shown by the time evolution of the system (Fig. |2J>c) . If the coupling p is far enough from this 
region, i.e., p < 0.25 or p > 0.4, the absorbent region loses its influence on the global dynamics and less structured 
and more random patterns than before are obtained (Fig. |2t--d). 

If the LMC complexity is quantified as function of p, our intuition is confirmed. The method proposed in Ref. pH | 
to calculate C is now adapted to the case of two-dimensional patterns. We let the system evolve until the stationary 
regime is attained. The variable xf in each cell is successively transformed in a binary sequence (0 if xf < 0.5 and 
1 if x™ > 0.5) when i covers the lattice, i — 1, . . . , L, and when n is consecutively increased. This binary string is 
analyzed in blocks of n bits, where n can be considered the scale of observation. The accessible states to the system 
among the 2™° possible states is found as well as their probabilities. Then, the magnitudes H, D and C are directly 
calculated. Figure El shows the result for the case of n a = 10. Let us observe that the highest C is reached when the 
dynamics displays spatio-temporal intermittency, that is, the most complex patterns are obtained for those values of 
p that are located on the borders of the absorbent region 0.258 < p < 0.335. Thus the plot of C versus p shows two 
tight peaks around the values p = 0.256 and p = 0.34 (Fig. |3J). 

If the detection of complexity in the two-dimensional case requires to identify some sharp change when comparing 
different patterns, those regions in the parameter space where an abrupt transition happens should be explored in 
order to obtain the most complex patterns. Smoothness seems not to be at the origin of complexity. As well as a 
selected few distinct molecules among all the possible are in the basis of life [28| , discreteness and its spiky appearance 
could indicate the way towards complexity. Let us recall that the distributions with the highest LMC complexity 
are just those distributions with a spiky- like appearance I n this line, the striking result here exposed confirms 
the capability of the LMC complexity for signaling a transition to complex behavior when regarding two-dimensional 
patterns. 
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FIG. 1: (x) versus p. The /3-statistics (or BW density) for each p is the rate between the number of black and white cells 
depicted by the system in the two-dimensional representation of its after-transient time evolution. (Computations have been 
performed with Ap = 0.01 for a lattice of 10000 sites after a transient of 5000 iterations and a running of other 5000 iterations). 
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FIG. 2: Digitalized plot of the one-dimensional coupled map lattice (axe OX) evolving in time (axe OY) according to Eq. (|2J: 
if Xi > 0.5 the (i, n)-cell is put in white color and if xf < 0.5 the (i, n)-cell is put in black color. The discrete time n is reset 
to zero after the transitory. (Lattices of 300 x 300 sites, i.e., < i < 300 and < n < 300). 



6 




FIG. 3: (□) C versus p. Observe the peaks of the LMC complexity located just on the borders of the absorbent region 
0.258 < p < 0.335, where /3 = (x). (Computations have been performed with Ap = 0.01 for a lattice of 10000 sites after a 
transient of 5000 iterations and a running of other 5000 iterations) . 



